Observation of Phase Separation in a Strongly-Interacting Imbalanced Fermi Gas 
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We have observed phase separation between the superfluid and the normal component in a strongly 
interacting Fermi gas with imbalanced spin populations. The in situ distribution of the density 
difference between two trapped spin components is obtained using phase- contrast imaging and 3D 
image reconstruction. A shell structure is clearly identified where the superfluid region of equal 
densities is surrounded by a normal gas of unequal densities. The phase transition induces a dramatic 
change in the density profiles as excess fermions start being expelled from the superfluid. 



PACS numbers: 03.75.Ss, 03.75.Hh, 05.70.Fh 

Cooper pairing is the underlying mechanism for the 
Bardeen- Cooper- Schrieffer (BCS) superfluid state of an 
equal mixture of two fermionic gases. An interesting sit- 
uation arises when the two components have unequal 
populations. Does the imbalance quench superfluidity, 
does it lead to phase separation between a balanced and 
imbalanced region, or does it give rise to new forms of 
superfluidity? A search for exotic superfluid states is 
promising in imbalanced mixtures, since the imbalance 
destabilizes BCS-type 5- wave pairing which is usually the 
strongest pairing mechanism 1, 2, 3, 4]. Recently, this 
problem has been experimentally addressed in ultracold 
atomic Fermi clouds with controlled population imbal- 
ances [3, 0, • Superfluidity was observed in a strongly 
interacting regime with a broad range of imbalances and 
the Pauli or Clogston limit of superfluidity 8] was char- 
acterized 

The phase separation scenario suggests that unpaired 
fermions are spatially separated from a BCS-superfluid of 
equal densities due to the pairing gap in the superfluid 
region |^, E|- our previous experiments 0, 0, 
we observed a strong central depletion in the difference 
profiles of expanding clouds indicating that excess atoms 
are expelled from the superfluid region. Ref. reports 
depletion of excess atoms at the trap center. It has been 
clarified 1^, J^] that none of these experiments answered 
the question: whether the densities of the two spin com- 
ponents are equal in the superfluid region and whether 
phase separation or rather distortions of the cloud due to 
interactions have occurred. 

Here we report the direct observation of phase sepa- 
ration between the superfluid and the normal region in 
a strongly interacting Fermi gas with imbalanced spin 
populations. The density difference between the two spin 
components is directly measured in situ using a special 
phase-contrast imaging technique and 3D image recon- 
struction. We clearly identify a shell structure in an im- 
balanced Fermi gas where the superfluid region of equal 
densities is surrounded by a normal gas of unequal den- 
sities. This phase separation is observed throughout the 
strongly interacting regime near a Feshbach resonance. 




FIG. 1: Phase-contrast imaging of the density difference of 
two spin states, (a) The probe beam is tuned to the red for 
the |1) \e) transition and to the blue for the |2) \e) 
transition. The resulting optical signal in the phase-contrast 
image is proportional to the density difference = ni — n2, 
where ni and n2 are the densities of the states |1) and |2), 
respectively, (b) Phase-contrast images of trapped atomic 
clouds in state |1) (left) and state |2) (right), and of an equal 
mixture of the two states (middle). 



Furthermore, we characterize the normal-to-superfluid 
phase transition of an imbalanced Fermi mixture using 
in situ phase-contrast imaging. The onset of superfluid- 
ity induces a dramatic change in the density profiles as 
excess fermions are expelled from the superfluid. 

A degenerate Fermi gas of spin-polarized ^Li atoms 
was prepared in an optical trap after laser cooling and 
sympathetic cooling with sodium atoms [illii. The 
population imbalance 5 of the two lowest hyperfine states 
|1) and 1 2) was adjusted with a radio- frequency sweep j^l- 
Here, 8 {Ni - N2) / {Ni ^ N2) , where TVi and 7V2 are the 
atom numbers in |1) and |2), respectively. Interactions 
between these two states were strongly enhanced near a 
broad Feshbach resonance at Bq = 834 G. The final evap- 
orative cooling was performed at 5 = 780 G by lowering 
the trap depth. Subsequently, the interaction strength 
was adiabatically changed to a target value by adjusting 
the value of the magnetic-bias field B with a ramp speed 
of < 0.4 G/ms. For typical conditions, the total atom 
number was Nt = Ni -\- N2 ~ 1 x 10^ and the radial 
(axial) trap frequency was fr = 130 Hz (/^ = 23 Hz). 

The condensate fraction in the imbalanced Fermi 



2 



] iitiiii 

5=5% 19% 27% 35% 41% 53% 62% 71% 75% 82% 88% 



FIG. 2: In situ direct imaging of trapped Fermi gases with 
population imbalance. The integrated 2D distributions of 
the density difference hd{x,z) = j nd(r) dy were measured 
using phase- contrast imaging. The data were obtained at 
B — 834 G for total atom number A^t 1 x 10^ and various 
population imbalance 5. For 5 < 75%, a distinctive core was 
observed showing the shell structure of the cloud. The field of 
view for each image is 160 /xm x 800 /xm. The three left-most 
images are displayed with different contrast levels for clarity. 
The image with 5 = 5(3)% was taken for Nt ^ 1.7 x lO'^. 



mixture was determined via the rapid transfer tech- 
nique |l^ [l^- Immediately after turning off the trap, 
the magnetic field was quickly ramped to B = 690 G 
{l/kpa ^ 2.6, where kp is defined as the Fermi momen- 
tum of a non-interacting equal mixture with the same 
total atom number and a is the scattering length) in ap- 
proximately 130 /is. The density profile of the expanding 
minority cloud was fit by a Gaussian for normal com- 
ponents (thermal molecules and unpaired atoms) and a 
Thomas- Fermi (TF) profile for the condensate 18]. 

The density difference between the two components 
was directly measured using a phase-contrast imaging 
technique described in Fig. ^ Iii this imaging scheme, 
the signs of the phase shifts due to the presence of atoms 
in each state are opposite so that the resulting phase sig- 
nal is proportional to the density difference of the two 
states if the probe frequency is adjusted properly [l^ . 
This technique allows us to directly image the in situ 
distribution of the density difference nd{r) between the 
two coniponents and avoid the shortcoming of previous 
studies 13, 0, where two images were subtracted from 
each other. 

For a partially superfluid imbalanced mixture, a shell 
structure was observed in the in situ phase-contrast im- 
ages (Fig. 12). Since the image shows the column density 
difference (the 3D density difference integrated along the 
^/-direction of the imaging beam), the observed depletion 
in the center indicates a 3D shell structure with even 
stronger depletion in the central region. The size of this 
inner core decreases for increasing imbalance and the core 
shows a distinctive boundary until it disappears for large 
imbalance. We observe this shell structure even for very 
small imbalances, down to 5(3)%, which excludes a ho- 
mogeneous superfluid state at this low imbalance, con- 
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FIG. 3: Reconstruction of 3D distributions from their inte- 
grated 2D distributions, (a) An integrated 2D distribution 
fid with 6 = 58% at B = 834 G. (b) A less noisy distribution 
was obtained by averaging four quadrants with respect to the 
central dashed lines in (a), (d) The ID profile obtained by 
integrating the averaged distribution along z-direction shows 
the 'flat-top' feature, (c) The 2D-cut, nd{x,y = 0,z), of the 
3D distribution nd{f) was reconstructed by applying the in- 
verse Abel transformation to the averaged distribution, as- 
suming only cylindrical symmetry along the 2;-direction. (e) 
The radial profile of the central section of the reconstructed 
3D distribution in xy-plane. The dashed lines in (d) and (e) 
are fits to the profiles' wings using a Thomas- Fermi (TF) dis- 
tribution. 



trary to the conclusions in Ref. |3| . 

The reconstructed 3D profile of the density difference 
shows that the two components in the core region have 
equal densities. We reconstruct 3D profiles from the 
2D distributions hd{x^z) of the column density differ- 
ence using the inverse Abel transformation (Fig.|3|) H^- 
The only assumption employed in this process is that 
of cylindrical symmetry of our trap along the axial z- 
direction. The two transverse trap frequencies are equal 
to better than 2% 15]. The reconstruction does not de- 
pend on the validity of the local density approximation 
(LDA) or a harmonic approximation for the trapping po- 
tential 0,0,113. ^ ID profile obtained by integrating 
hd along the axial z-direction (Fig.|3fd)) shows a flat top 
distribution, which is the expected outcome for a shell 
structure with an empty inner region in a harmonic trap 
and assuming LDA. 

The presence of a core region with equal densities for 
the two components was correlated with the presence of 
a pair condensate. The density difference at the cen- 
ter Udo along with the condensate fraction is shown as 
a function of the population imbalance S in Fig. ^ As 
shown, there is a critical imbalance 6c where superfluid- 
ity breaks down due to large imbalance [EH- In the 
superfluid region, i.e., S < Sc, rido vanishes, and for 
S > Sc^ rido rapidly increases with a sudden jump around 
S ^ 5c. We observe a similar behavior throughout the 
strongly interacting regime near the Feshbach resonance, 
—0.4 < l/kpci < 0.6. This observation clearly demon- 
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FIG. 4: (Color online) Phase separation in a strongly interact- 
ing Fermi gas. Normalized central density difference 77 (black 
circle) and condensate fraction (red triangle) as functions of 
imbalance S for various interaction strengths. 77 = Udo/no 
where Udo was measured as the average over the central region 
of 7 /xm x40 /im in a 2D-cut of the 3D distribution (Fig.^Ic)) 
and no = 3.3 X 10^^ cm~^ is the calculated central density of 
a fully polarized Fermi gas with A^t = 1 x 10^. The conden- 
sate fraction was determined from the bimodal distribution 
of the minority component after a magnetic field sweep and 
expansion. The solid line is a fit for the condensate fraction 
to a threshold function oc {l — \6/6c\^). The critical imbalance 
6c indicated by the vertical dashed lines were 96%, 77%, and 
51% for B = 780 G, B = 834 G and B = 884 G, respectively. 



strates that for this range of interactions a paired super- 
fluid is spatially separated from a normal component of 
unequal densities. 

The shell structure is characterized by the radius of the 
majority component, the radial position Rp of the peak 
in nd{r)^ the size Rc of the region where Ud is depleted, 
and the "visibility" of the core region a = {nd{Rp) — 
^do)/(^d(^p) + ^do) (Fig. [5|). The sudden drop of a 
around S ^ Sc has the same origin as the observed sudden 
jump of UdO' the breakdown of superfluidity (compare 
Fig.^. The comparison between Rp and Rc shows that 
the boundary layer between the superfluid and normal 
region is rather thin. It has been suggested that the de- 
tailed shape of profiles in the intermediate region could 
be used to identify exotic states such as the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state HHHIl^. This wih 
be a subject of future research. 

The superfluid requires equal central densities in the 
strongly interacting regime at our coldest temperatures. 
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FIG. 5: (Color online) Characterization of the shell structure, 
(a) Reconstructed 3D radial profiles at B = 834 G. (b) Radius 
of the majority component R (black circle), peak position of 
the density difference Rp (red circle) , and radius of the empty 
core Rc (blue triangle). R was determined from the profiles' 
wings using a fit to a zero-temperature TF distribution. The 
solid line indicates the TF radius of the majority component 
in an ideal non-interacting case. R ai 6 = (open circle) 
was measured from an image taken with a probe frequency 
preferentially tuned to one component. Rc is defined as the 
position of the half-peak value in the empty core region for 
6 < 6c. (c) Visibility of the core region is defined as a = 
{nd{Rp) - rido) / {nd{Rp) + Udo). 



A normal imbalanced Fermi mixture will have unequal 
densities. Thus, one should expect that a visible change 
in the density difference occurs as the temperature is 
lowered across the normal-to-superfluid phase transi- 
tion [3, |9|] . In situ phase-contrast images of a cloud at 
various temperatures are shown in Fig. El The tempera- 
ture T of the cloud is controlled with the final value of the 
trap depth in the evaporation process. The shell struc- 
ture appears and becomes prominent when T decreases 
below a certain critical value. This shell structure gives 
rise to the bimodal density profile of the minority com- 
ponent that we observed after expansion from the trap 
in our recent work . Here we show via in situ measure- 
ments that the onset of superfluidity is accompanied by 
a pronounced change in the spatial density difference. 

The phase transition is characterized in Fig. [3 As T 
is lowered, Udo gradually decreases from its plateau value 
and the condensate fraction starts to increase. From the 
point of condensation (condensate fraction > 1%) and 
deformation of the minority clouds (x^ in Fig.[7Jb)), we 
determine the critical temperature Tc = 0.13(2) Tp for 
the imbalance of ^ = 56(3)%. Tp = 1.7 /iK is the Fermi 
temperature of a non-interacting equal mixture with the 
same total atom number. The rise in the drop in rido 
and the onset of condensation are all observed at about 
the same temperature. Better statistics are needed to 
address the question: whether some weak expulsion of 
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FIG. 6: Emergence of phase separation in an imbalanced 
Fermi gas. The temperature of the cloud was controhed by 
varying the final value of the trap depth Uf in the evaporation 
process. Phase- contrast images were taken after adiabatically 
ramping the trap depth up to ks x 3.7 /iK {fr = 192 Hz). 
The whole evaporation and imaging process were performed 
at 5 = 834 G{Nt^ 1.7x10^, 6^ 56%). The field of view for 
each image is 160 /xm x 940 /xm. The vertical and horizontal 
scale of the images differ by a factor of 1.5. 
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FIG. 7: (Color online) Phase transition in an imbalanced 
Fermi gas. The phase transition shown in Fig. |S] was charac- 
terized by measuring (a) population imbalance S, (b) temper- 
ature T (black circle), for fitting the minority cloud with a 
finite temperature Fermi-Dirac distribution (red triangle), (c) 
central density difference Udo (black circle), and condensate 
fraction (red triangle). The solid line is a guide to the eye 
for x^- When Uf < /cs x 1 /iK, S decreased mainly due to 
loss of the majority component. The number of the minority 
component was constant (A^2 ~ 3.7 x 10^). T was determined 
from the non-interacting outer region of the majority cloud 
after 10 ms of ballistic expansion 7]. T and n^o were averaged 
values of three independent measurements. The open (solid) 
arrow in (c) indicates the position for Tc (T*). See the text 
for the definitions of Tc and T*. 



majority atoms from the center occurs already slightly 
above Tc. 

Below a certain temperature T* , Udo reaches zero while 
the condensate fraction is still increasing, implying that 
the superfluid region of equal densities continues to ex- 
pand spatially with decreasing T. Full phase separa- 
tion does not occur until this temperature T* < Tc 
is reached. We interpret the state between T* and Tc 
as a superfluid of pairs coexisting with polarized quasi- 
particle excitations. This is expected, since at finite tem- 
peratures the BCS state of an equal mixture can accom- 
modate excess atoms as fermionic quasi-particle excita- 
tions 0,12^ . There is a finite energy given by the pairing 
gap A(T) for those quasi-particles to exist in the super- 
fluid. The assumption that excess atoms should have 
thermal energy ksT > A(T) to penetrate the super- 
fluid region suggests the relation between T* and A(T) 
to be ksT^ ^ A(T*). From our experimental results, 
T* ^ 0.09 Tf and A(T*) ^ h x 3.3 kHz. 

In conclusion, we have observed phase separation of 
the superfluid and the normal component in a strongly- 
interacting imbalanced Fermi gas. The shell structure 
consisting of a superfluid core of equal densities sur- 
rounded by a normal component of unequal densities was 
clearly identified using in situ phase-contrast imaging and 
3D image reconstruction. The phase-contrast imaging 
technique combined with the 3D reconstruction process 
provides a new method to measure the in situ density 
distribution, allowing direct comparison with theoretical 
predictions. 
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